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Self-Ignition Temperatures of Materials From Kinetic- 
Reaction Data 

Daniel Gross and A. F. Robertson 

Results of experimental determinations of the kinetic constants of the self-heating 
reaction are presented for wood fiberboard, cotton linters, sugar pine, cork, crepe rubber, 
GRS rubber, natural, synthetic, and blended foam rubber (with and without additive), 
various oils (raw linseed, cottonseed, rapeseed, sperm, olive, castor, and neatsfoot) applied 
to cotton gauze in a ratio of 1 part of oil to 6 parts of cotton by weight, ammonium per- 
chlorate, and nitrocellulose plastic. Under the assumption that self-heating follows a 
first-order reaction, these constants were used to calculate the critical radii of spherical 
piles for each of four surface temperatures likely to be experienced in long-period storage. 
Calculated self-ignition temperatures of piles of ^-inch-diameter to 22-inch-diameter spheres 
of wood fiberboard and %-inch-diameter to 2-inch-diameter spheres of cotton linters were 
in reasonable agreement with previous measurements by N. D. Mitchell (National Fire 
Protection Association Quarterly 45: 162, 1951). 

1. Introduction 

Previous workers studying the self-ignition be- 
havior of materials have experienced difficulties in 
the determination and specification of a unique 
ignition temperature. Because of the critical heat 
balance that exists during the ignition process, the 
measured "ignition temperature" is dependent upon 
the characteristics of the testing apparatus, the 
degree of control of the ambient temperature condi- 
tions, the time of exposure, and other experimental 
factors. It is highly dependent also upon the size 
of the specimen under test. This is due to the fact 
that self-ignition of a material is only possible when 
the rate of heat generation within the material ex- 
ceeds the rate at which heat is lost by radiation, 
convection, and conduction to the surroundings. 
Since the heat generated is a function of the volume 
of the material, while the heat lost is a function of 
its surface area, the measured ignition temperature 
should decrease as the linear dimensions increase. 

The incidence of fires resulting from self -ignition 
makes desirable an understanding of the relationship 
between the self-heating reaction and the physical 
parameters involved. The present paper proposes 
and illustrates the manner in which kinetic-reaction 
data obtained under conditions of negligible heat 
loss may be used to predict critical surface and 
ambient temperatures as a function of size and shape 
for a variety of materials. 



2. Analysis and Nomenclature 

In the analysis, the following symbols and units 
will be used: 



A heat-generation coefficient 

B radius or half-thickness 

B c critical radius of sphere 

c specific heat 

e Napierian base 

E activation energy 

h surface heat-transfer co- 
efficient 

k thermal conductivity 

i)i coefficient in eq (1) 

R gas constant =0.0019872 



cal sec _1 cm~ 3 

cm 

cm 

cal g-i° K-i 



kcal mole -1 

cal sec-icm- 20 K" 1 

cal sec -1 cm~ 10 K -1 

kcalmole-^K 11 "" 



t 


i [me 


sec 


T 


temperature within mate- 
rial 


°K 


T A 


self-ignition (critical am- 
bient) temperature 


°K 


T B 


critical surface tempera- 
ture 


°K 


To 


critical center temperature 


°K 


T m 


mean temperature = ( To + 
Tb)I2 


°K 


X 


space coordinate 


cm 


p 


density 


gen 



Analyses of self-heating have been made under the 
assumption that heat is generated within the body 
according to a first-order reaction law. Although 
this may be an oversimplification for complex mate- 
rials and over extended temperature ranges, the 
assumption appears to be applicable for many mate- 
rials in the relatively low-temperature ranges where 
self-heating becomes important and the ignition 
process is initiated. Provided the rate of the self- 
heating process within the material is not limited by 
diffusion of the reactants and/or products, or by 
depletion of the fuel supply, a steady-state condi- 
tion can exist when the heat generated within the 
material is equal to the heat lost to the surroundings. 
The particular condition of interest is the one that 
produces the maximum self-heating which is still 
steady state. Either a slightly higher ambient tem- 
perature or a slightly larger bulk size results in a 
greater quantity of heat being generated than can 
be dissipated, and ignition results. Since this steady- 
state condition is critical with respect to ignition, the 
corresponding surface temperature and size of mate- 
rial are referred to as the critical surface temperature 
and the critical size. The ambient temperature cor- 
responding to given critical steady-state conditions 
with given heat losses is considered the critical ambi- 
ent or self -ignition temperature. There exists for 
each material and for each set of heat-loss conditions 
a series of corresponding self-ignition temperatures 
and critical sizes. 

Consider a homogeneous material of half-thickness 
B generating heat in accordance with the first-order 
reaction law Aexp~E/RT and losing heat from its 
surface. The steady-state heat-conduction equation 
which applies for the condition when the heat lost 



413 



through the surface is equal to the heat generated 
within the material may be written : 

The critical condition of ignition in the theory of 
thermal explosions originally proposed by Frank- 
Kamenetzky [1] 1 requires the solution of this non- 
linear equation. Chambre [2] presented analytic 
solutions in terms of known functions for the sphere 
(m=2) and cylinder (m=l) to supplement that for 
the semi-infinite slab (m=0). Recent analyses by 
Genensky [3] for a sphere and by Thomas [4] for 
all three geometries have been made employing the 
assumption common to the previous solutions, 
namely, that temperature differences within the ma- 
terial are small in comparison with the temperature. 
Enig, Shanks, and Southworth [5] have provided 
solutions in the form of convenient tables of the 
center and surface temperatures for sphere, cylinder, 
and semi-infinite slab geometries without this restric- 
tive assumption. From the tables it is useful to 
note that for a given surface temperature T B , the 
critical radius for the cylinder and the critical half- 
thickness for the semi-infinite slab are given very 
closely by 0.775 B c and 0.514 B c , respectively, where 
B c is the critical radius for the sphere. Assuming 
that the heat generated within the material at a 
mean temperature T m is lost from its surface accord- 
ing to Newton's law of cooling, the following relation 
may be used for calculating the critical ambient 
temperature for a sphere : 



T A =T B 



AB C - 



3h 



E 



(2) 



where T m is a mean temperature, here assumed as 
To+T B l2. 

Since the analysis by Thomas employs Newtonian 
surface cooling as a boundary condition, the ambient 
temperature is obtained directly. From the graphi- 
cal relations given, it is again useful to note that for 
a given ambient temperature T A , the critical radius 
for the cylinder and the critical half-thickness for the 
semi-infinite slab are given very closely by 0.782 B c 
and 0.520 B c respectively. 

The experimental method used in this study for 
determination of the constants of the first-order 
reaction involved self-heating conditions which 
differed from those of the self-heating process 
described. Instead of a fixed ambient temperature, 
the temperature within the adiabatic test furnace 
surrounding the specimen was controlled to maintain 
negligible loss of heat. Under these conditions, the 
specimen temperature remained uniform throughout 
its mass and the heat generated within the body 
increased its temperature according to the relation: 



dT -— 
pc-Tr=Ae RT > 
at 



(3) 



By transposing terms, taking logarithms and 
plotting In dTjdt versus 1/T, the resultant line has a 
slope of— (E/R)&nd intercepts the In dTjdt axis at In 
A/pc. For each material the determination of the 
values of A and E, as well as of its more common 
thermal properties, permits the solution of the 
general self-heating problem including the case 
where size and temperature become critical and 
ignition results. 

3. Experimental Apparatus 

In the determination of the kinetic constants, use 
was made of an improved version, figures 1 and 2, of 
an adiabatic furnace previously described [6]. 
Although the basic operation remains the same, 
refinements have been made and are here sum- 
marized : 

Addition of guard heaters and associated control 
system. In order to reduce heat losses and to main- 
tain minimum temperature gradients within the 
furnace enclosure, the Dewar flask has been enclosed 
within a close-fitting cylindrical shell. The shell 
consists of two concentric stainless steel cylinders 
containing two electric heating elements and an 
insulating fill. The shell contains the bottom and 
center guard heating elements while a top-guard 
heating element has been placed within the Dewar 
flask plug composed of several layers of asbestos 
board. A thyratron control system as shown in 
figure 3 permits adjustment of the guard-heater 
cycle. 

Use oj multiple-junction thermocouple. The sens- 
ing element for the servo controller consists of eight 
No. 30 gage (B&S) copper-constantan thermo- 
j unctions arranged in series to measure the mean 
temperature-difference between the specimen and 
its ambient medium. 




1 Figures in brackets indicate"the literature reference at the end of this paper. 



Figure 1. Adiabatic furnace with control and recording 
equipment. 
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Figure 2. Adiabatic furnace, showing Dewar flask and speci- 
men container removed from insulating shell. 



The changes and additions described permit im- 
proved performance of the system. The controller 
can sense temperature differences of the order of 
0.01° C while the rate of temperature drift as ob- 
served with an inert specimen has been reduced to 
approximately 0.25° C/hr over the temperature 
range 30° to 300° C. 

4. Results 

4.1. Determination of A and E 

Table 1 lists some thermal and kinetic constants 
for a number of materials. Some of the values of k 
and c were determined during this investigation, 
while others were obtained from handbook sources. 
The values of A and E were determined over the 
indicated temperature ranges from experimental 
temperature- versus-time data and the use of eq (3). 
For some of these materials, there is reasonable agree- 
ment with published values obtained by other 
methods. These values were determined over lim- 
ited temperature ranges before loss of weight became 
appreciable. An analysis that includes the effect of 
weight loss may be applied for the higher tempera- 
ture ranges. Analysis of several cellulosic materials 
has actually shown that when weight changes are 
taken into consideration, the values of the kinetic 
constants were not much different from those 
obtained at the lower temperature ranges. 



SWITCHES 

51 DPST MAIN POWER 

52 SPST RESET 

53 DPST CONSTANT 
TEMPERATURE CONTROL 

54 DPST MAIN HEATER 

55 SPST VOLTMETER 

56 SPST LIMIT (MERCURY) 



/ S4 / 




R I, R9 TIME DELAY 
R2,R3DPDT LATCHING 
R4 SPST CONTROL 
R5 SPST VOLTMETER 
R6 SPST LIMIT (MERCURY) 
R7, R8 SPST THERMOSTATIC 
CONTROL 



SERVO 
MOTOR 



Figure 3. Circuit diagram of furnace-control jy stem. 
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Table 1. Physical and kinetic properties 



Material 



Temperature 
range 



Property 



k 


E 


cnlsec- 1 cm- ] K~ l 


kcal mole~ l 


0. 00012 


25.7 


. 00010 


34.5 


. 00024 


20.3 


. 00010 


19.1 


. 000282 


23.4 


. 000299 


23.6 


. 000096 


27.6 


. 000100 


26.2 


. 000098 


29.0 


. 000095 


28.5 


.00011 


21.0 


. 00011 


24.1 


. 00011 


20.1 


. 00011 


20.1 


. 00011 


18.7 


. 00011 


27.8 


. 00011 


24.7 




41.2 
42.0 


. 00051 



Wood fiberboard 

Cotton linters 

Sugar pine 

Cork 

Crepe rubber 

GRS 

Foam rubber, 100% natural 

Foam rubber, 100% synthetic 

Foam rubber, blended, with additive 

Foam rubber, blended, without additive. 

Raw linseed oil * 

Cottonseed oil a 

Rapeseed oil a 

Sperm oil » 

Olive oil a 

Castor oil a 

Neatsfoot oil a 

Ammonium perchlorate 

Nitrocellulose plastic 



°C 
80-225 
180-260 
140-310 
170-225 
210-310 
220-330 
100-180 
115-210 
110-190 
110-200 
90-120 
75-120 
90-140 
120-160 
115-160 
140-160 
140-170 
205-270 
135-170 



g cm~* 
0.25 
.32 
.36 
.13 
.92 
1.10 
.108 
.129 
.118 
.104 
.309 
.338 
.316 
.313 
.308 
.314 
.316 
.91 
.82 



calg^°K 
0.33 
.32 
.32 
.48 
.5 
.5 
.5 
.5 
.5 
.5 
.34 
.34 
.34 
.34 
.34 
.34 
.34 
.35 
.34 



cal sec~ l cm~' 
1.97X108 
5.30X10" 
2.11X106 
2.50X105 
1.06X106 
1.11X107 
7. 48X10io 
4.96X109 
4.25X10" 
1.10X10H 
3.50X109 
2.42X10" 
6. 13X108 
1.56X107 
2.41X107 
8.33X10" 
3. 82X1010 
1.22X1016 
2.43X1019 



a Applied to cotton gauze in ratio of 1 part of oil to 6 parts of cotton by weight. All property values refer to oil-cotton combination. 



4.2. Determination of Critical Size and Temperature 

The critical surface temperatures of spherical piles 
of wood fiberboard and cotton linters were computed, 
using the measured values of the kinetic and thermal 
constants, and the tables given in reference [5]. 
These computations were performed over the radius 
range corresponding to the specimen sizes reported 
by Mitchell [7]. Equation (2) was then used to cal- 
culate the corresponding critical ambient tempera- 
tures. The value of the surface heat-transfer 
coefficient was assumed as 4.4XKH cal/sec cm 2 °K. 
Using this value for the surface heat-transfer coeffi- 
cient, the critical ambient temperature was also 
determined from the relations given in reference [4] . 
The results of these calculations are shown in table 2 
for two materials for which published experimental 
data [7] are available. With the exception of the 
smallest size specimens, for which the assumption of 
an arithmetic mean temperature for eq (2) introduces 

Table 2. Comparison between measured and calculated 
self-ignition temperatures, Ta 



Size 


Wood fiberboard 


Cotton linters 


Thick- 
ness 
2B e 


Half- 
thick- 
ness 
B e 


Self-ignition temperature 


Self-ignition temperature 


Meas- 
ured a 

Ref7 


Calculated b 


Meas- 
ured a 


Calculated b 


Ref 5 and 
eq2 


Ref 4 


Ref 7 


Ref 5 and 
eq2 


Ref 4 


in. 

H 
H 

l 

2 
4 
8 
12 
22 


cm 
.159 
.318 
.635 
1.27 
2.54 
5.08 
10.2 
15.2 
27.9 


°C 
317 
252 
212 
202 
180 
150 
130 
122 
109 


°C 

183 
218 
215 
200 
180 
159 
140 
129 
114 


°C 
260 
243 
223 
202 
180 
160 
140 
129 
114 


°C 
304 
270 
248 
222 
201 


°C 
257 
274 
267 
251 
233 


°C 
307 
291 
272 
253 
233 



a Specimens \i, H, and Y in. thick were cubes; specimens 1 in. and larger 
were octagonal prisms. 
*> Based on spheres of radius equal to half-thickness. 



considerable error in the calculated self -ignition 
temperature, the results are seen to be in reasonable 
agreement. 

Table 3 lists the critical radius B c of a sphere corre- 
sponding to several critical surface temperatures for 
each of the materials listed in table 1. Curves have 
been plotted in figure 4 showing this relationship for 
seven materials. In this temperature range (20° to 
100° C), of practical interest for ordinary storage, the 
critical size calculations were based upon the assump- 
tion of the applicability of the kinetic constants which 
were determined over higher temperature ranges. It 
was further assumed that the thermal properties re- 
mained constant over the whole temperature range. 
Critical dimensions for cylinders and slabs may [be 
estimated from the previously mentioned constants of 
0.77 bB c and 0.514B C , respectively. As h becomes very 
large, the self-ignition temperature is equal to the 

Table 3. Computed critical radius of a sphere, B c 



Material 


Critical surface temperature 


20° C 


48.9° C 


82.2° C 


100° C 


Wood fiberboard 


cm 
4.49X103 
4.25X106 
2.19X103 
1.43X103 
1.43X10* 
1.71X10* 

3.36X103 

4.00X103 

4.65X103 

4.53X103 
66 
101 
72 
143 
110 
1.16X103 

139 
7.69X10* 


cm 

680 
3.29X105 

501 

363 
2.58X103 
3.03X103 

438 

588 
542 

565 
14.3 
17.3 
17.0 
34 
29 

150 

37 

3.31X103 


cm 
114 
2.86X10* 
125 
99 
511 
591 

64 

95 

71 

78 
3.6 
3.3 
4.3 
8.5 
7.3 

69 
4.1 
167 


cm 

50 
9.43X103 

66 

55 
206 
280 

26 

41 

27 

34 
1.7 
1.5 
2.3 
4.5 
4.5 
9.0 
1.9 

43 




Cork. 




GRS rubber 


Foam rubber, 100% 
natural- -- - 


Foam rubber, 100% 
synthetic. .. - 


Foam rubber, blended, 
with additive 

Foam rubber, blended, 
without additive 

Raw linseed oil a 

Cottonseed oil a 




Olive oil a __ - -- - - 






Nitrocellulose plastic. .. 



a Applied to cotton gauze in ratio of 1 part of oil to 6 parts of cotton" by 
weight. All property values refer to oil-cotton combination. 
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20 40 60 80 100 

CRITICAL SURFACE TEMPERATURE ,°C 

Figure 4. Computed critical radii and critical-surface tem- 
peratures for spherical piles of various materials. 



critical surface temperature. For a finite h, the self- 
ignition temperature may be estimated from eq (2) or 
may be obtained directly from the relations presented 
in reference [4]. 

5. Conclusions 

Recent analyses of self -ignition have been made 
under the assumption that self -heating follows a 
first-order reaction law. Experimental determina- 
tions of the kinetic constants of this reaction have 
been made for a number of materials and these con- 
stants have been used to calculate self-ignition 
(critical ambient) temperatures as a function of size 
and shape. The calculated self-ignition temperatures 
are in fair agreement with published experimental 
data for two materials, wood fiberboard, and cotton 
linters. 



The experimental determinations and most of 
the computations were performed by Joseph J. 
Loftus to whom the authors are deeply appreciative. 
Thanks are due also to James S. Evans, Jr., who 
designed the thyratron control circuits and to 
William H. Bailey, who assembled the furnace and 
made the drawings. 

An important contribution to this field has been 
made by Kazuo Akita of the Fire Research Institute, 
Tokyo, in a paper to be published in J. Chem. Soc. 
Japan, Section of Indust. Chem. 
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